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SUMMARY 


The  purpose  of  the  theoretical  and  experimental  study  reported  herein 
was  to  investigate  the  flow  characteristics  of  circular  uniform  jets  imping¬ 
ing  normally  on  a  ground  plane  in  order  to  further  the  understanding  of 
the  aerodynamic  processes  associated  with  ground  particle  entrainment 
in  the  impinging  downwash  of  VTOL  aircraft.  Particular  emphasis  was 
placed  on  the  flow  properties  in  the  ground-plane  boundary  layer  where 
entrainment  occurs.  On  the  basis  of  the  investigation  that  was  conducted, 
the  following  general  conclusions  were  drawn. 


Experimental 


The  effect  of  the  ground  on  the  jet  exit  velocity  profile  for  a  circular 
uniform  jet  impinging  normally  on  the  ground  is  negligible  for  jet  exits 
more  than  1.  5  to  2.0  nozzle  diameters  above  the  ground  plane.  For 
jet  heights  of  more  than  2.  0  nozzle  diameters  above  the  ground  plane, 
changes  in  the  jet  flow  characteristics  with  height  are  primarily  associ¬ 
ated  with  the  effects  of  viscous  mixing. 


Three  distinct  boundary  layer  regions  occur  in  the  ground-plane  flow. 

A  laminar  boundary  layer  exists  for  0  ^.r/D  ^0.  8;  a  transition  region 
of  mixed  turbulent  and  laminar  boundary  layer  flow  occurs  for 
0.  8  £>  r/Di  2.0;  beyond  r /D  =  2.0  for  values  of  the  nozzle  height-to- 
diameter  ratio,  H/D,  less  than  or  equal  to  4,  the  ground  flow  is  of  the 
turbulent  radial  wall-jet  type.  Characteristics  and  extent  of  the  transi¬ 
tion  boundary  layer  region  vary  with  H/D.  For  the  larger  values  of  H/D, 
the  transition  to  turbulent  flow  appears  to  be  forced  by  interaction  of 
the  mixing  region  of  the  jet  with  the  boundary  layer.  The  start  of  the 
fully  developed  wall-jet  flow  occurs  at  values  of  r/D  which  decrease  as 
H/D  decreases. 
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Theoretical 

Currently  available  methods  of  predicting  the  inviscid  three-dimensional 
impinging  jet  flow  suffer  from  the  lack  of  an  accurate  method  of  predicting 
the  shape  and  position  of  the  free -stream  surface.  A  plausible  iterative 
technique  is  proposed  for  doing  this,  but  numerical  results  using  this 
technique  have  not  yet  been  obtained.  Calculations  of  ground-plane  pres¬ 
sure-distribution  based  on  a  jet  boundary  derived  by  Schach  (Reference  8) 
for  the  infinite  impinging  jet  and  using  the  vortex  sheet  formulation  for 
the  jet  developed  during  the  present  research  resulted  in  relatively  good 
agreement  with  experiment. 

An  approximate  three-dimensional  boundary  layer  analysis  accurately 
predicts  the  laminar  portion  of  the  boundary  layer,  using  experimentally 
derived  pressure  gradients,  but  fails  to  predict  the  experimental  transi¬ 
tion  point  at  all  H/D's.  This  failure  is  apparently  a  result  of  the  effect 
of  the  jet  mixing  region  interacting  with  the  boundary  layer  prior  to 
natural  transition  for  H/D-£  1.  0.  Glauert' s  wall-jet  theory  is  directly 
applicable  to  the  ground-plane  flow  for  the  r/D>  2.  0  for  all  H/D's 
tested.  A  semiempirical  method  was  developed  for  predicting  boundary 
layer  characteristics  in  the  transition  region.  This  method  in  conjunc¬ 
tion  with  the  Glauert  theory  allows  prediction  of  velocities  and  boundary 
layer  thickness  to  within  about  10  percent  of  experimental  values  in 
the  wall-jet  flow  regime  of  the  ground  flow. 
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CONCLUSIONS 


On  the  basis  of  the  results  obtained,  the  following  conclusions  can  be 

drawn  regarding  the  experimental  flow  properties  of  uniform,  circular, 

normal  impinging  jets,  and  regarding  the  theoretical  analyses  conducted: 

1.  The  effect  of  the  ground  on  the  jet- exit  velocity  profile  is  neg¬ 
ligible  for  jet  exits  more  than  one  and  one-half  to  two  nozzle 
diameters  above  the  ground  plane. 

2.  Changes  in  the  jet  flow  characteristics  with  jet  nozzle  height 
are  primarily  associated  with  the  effects  of  viscous  mixing  for 
nozzle  heights  greater  than  two  nozzle  diameters  above  the 
ground  plane. 

3.  For  values  of  H/D<2  the  pressure  distribution  on  the  ground 
under  a  viscous  impinging  jet  is  essentially  that  obtained  with 
an  inviscid  jet. 

4.  Transition  from  laminar  to  turbulent  flow  in  the  ground-plane 
boundary  layer  is  considerably  affected  by  interaction  of  the 
jet  mixing  region  with  the  ground  boundary  layer. 

5.  The  start  of  a  radial  wall-jet  flow  occurs  at  radial  distances 
from  the  ground-plane  stagnation  point  which  decreases  with 
decreasing  jet  nozzle  height.  For  small  nozzle-to-ground  dis¬ 
tances  (H/D<2),  this  decrease  is  primarily  associated  with 
the  reduction  in  mean  velocity  in  the  nozzle  as  H/D  decreases. 

6.  Approximate  boundary  layer  analyses  resulted  in  good  to  excel¬ 
lent  agreement  with  experimental  data,  except  in  the  small 
region  between  the  radius  where  transition  to  turbulent  flow 
occurs  and  the  radius  at  which  the  velocity  profiles  first 
approach  the  form  predicted  by  Glauert,  for  which  region 
there  is  no  applicable  theory. 
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RECOMMENDATIONS 


On  the  basis  of  the  experimental  results  obtained  and  the  analyses  con¬ 
ducted,  the  following  recommendations  are  made: 

1.  An  experimental  investigation  of  the  jet  flow  characteristics 
and  boundary  layer  flow  characteristics  should  be  made  of  a 
normally  impinging  nonuniform  jet  representative  of  the  down- 
wash  of  helicopter  rotors  and  ducted  fans. 

2.  Theoretical  calculations  of  the  flow  field  of  the  inviscid  uniform 
impinging  jet  should  be  made  for  the  range  of  jet  nozzle  height - 
to-diameter  ratios  of  practical  interest. 

3.  Experimental  measurements  should  be  made  of  the  aerodynamic 
forces  on  objects  in  flows  representative  of  the  ground  flows 
generated  by  jets,  ducted  fans,  and  rotors. 
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I  INTRODUCTION 


Of  the  many  VTOL/STOL  operational  problems,  one  of  the  more  serious 
is  that  of  particle  entrainment  in  jet  or  rotor  downwash  during  near- 
ground  hovering  or  slow -speed  flight.  Such  entrainment  can  lead  to 
loss  of  visibility  for  the  pilot,  loss  of  concealment  for  the  aircraft,  and 
damage  to  powerplant  and  structure.  A  considerable  body  of  literature 
has  grown  over  the  past  few  years  concerning  various  aspects  of  this 
problem.  References  1,  2,  and  3  are  typical  of  the  experimental 
studies  of  impinging  downwash  from  jets,  ducted  fans  and  rotors;  data 
have  been  obtained  on  both  particle  entrainment  and  the  associated 
flow  field. 

Vidal,  in  Reference  4,  reviewed  the  impinging -downwash,  particle- 
entrainment  problem  from  the  point  of  view  of  available  experimental 
data  and  theory  in  order  to  determine  mechanisms  for  entrainment.  A 
drag  criterion  for  particle  entrainment  has  been  proposed  (Reference  1), 
based  on  the  drag  required  to  start  a  particle  moving  while  in  contact 
with  the  ground.  Particles  are  entrained  in  the  flow  by  bouncing  off 
projections  on  the  ground.  Vidal  proposed,  in  addition,  a  lift  criterion 
for  entrainment,  based  on  the  lift  required  to  pick  up  a  particle  clear 
of  the  ground.  The  lift  on  the  particle  is  dependent  on  the  velocity 
gradient  in  the  ground  boundary  layer,  as  well  as  on  the  mean  dynamic 
pressure  of  the  flow  at  the  particle.  On  the  basis  of  the  drag  and  lift 
criteria  for  entrainment,  and  using  certain  simplifying  assumptions, 
Vidal  showed  qualitative  agreement  between  calculated  areas  of  particle 
entrainment  under  an  impinging  downwash  and  areas  estimated  from 
available  experimental  data.  The  relative  importance  of  the  two  cri¬ 
teria  depends  primarily  on  the  size  of  the  particle  compared  to  the 
boundary  layer  thickness. 
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A  quantitative  understanding  of  the  entrainment  mechanism  may  well 
lead  to  effective  means  of  preventing  or  at  least  reducing  the  severity 
of  the  particle  entrainment  problem  in  impinging  downwash  flows.  Inas¬ 
much  as  entrainment  occurs  in  the  flow  along  the  ground  and  as  the 
boundary  layer  appears  to  play  an  important,  if  not  dominant,  part  in 
this  process,  the  need  for  accurate  experimental  data  for  the  ground- 
flow  boundary  layer  is  apparent.  Although  there  is  a  considerable  body 
of  experimental  data  for  both  impinging  jets  and  rotor  downwash  flows, 
to  the  authors'  knowledge  there  is  none  in  which  a  detailed  study  of  the 
boundary  layer  in  the  flow  along  the  ground  has  been  made.  The  experi¬ 
mental  and  theoretical  results  obtained  from  this  and  further  research 
should  provide  the  basis  for  the  development  of  an  accurate  analytical 
technique  for  predicting  areas  under  the  impinging  downwash  of  VTOL- 
type  aircraft  where  ground  particles  can  be  entrained.  It  is  believed 
that  the  development  of  effective  and  practical  methods  of  preventing 
or  at  least  alleviating  the  serious  consequences  of  particle  entrain¬ 
ment  would  be  materially  advanced  thereby. 

The  research  described  in  this  paper  was  an  experimental  and  analytical 
program  to  study  the  flow  processes  in  a  uniform  impinging  jet.  Ana¬ 
lytical  methods  for  predicting  the  flow  in  such  a  jet,  with  particular 
emphasis  on  the  boundary  layer  flow,  were  the  primary  aim.  Experi¬ 
mental  data  on  the  flow  in  the  uniform  impinging  jet  and  in  the  ground 
boundary  layer  are  presented,  as  are  the  results  of  the  theoretical 
studies  obtained  thus  far. 
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II  EXPERIMENTAL  PROGRAM 


Experimental'  Setup 

The  CAL-Air  Force  One-Foot  High-Speed  Wind  Tunnel  was  used  during 
the  experimental  research.  The  transonic -super sonic  test  section  was 
completely  removed.  A  12 -inch-diameter  molded  Fiberglas  nozzle, 
designed  to  produce  a  uniform  and  parallel  flow  at  the  exit,  was  mounted 
to  an  adapter  plate.  The  plate  was  attached  to  the  test-section  fa£.e  of  the 
wind  tunnel  settling  chamber.  The  ground  board  was  8  feet  x  8  feet 
square  and  could  be  set  at  varying  distances  from  the  jet  nozzle.  Fig¬ 
ures  2  and  3  show  the  experimental  setup.  A  pattern  of  static  pressure 
taps  was  located  in  the  ground  board  (Figure  4). 

Velocities  and  static  pressures  in  the  jet  were  measured  with  a  Pitot- 
static  probe  mounted  on  a  traverse  mechanism.  Ground  flow  measure¬ 
ments  were  obtained  from  static  and  total  head  probes  mounted  in  a 
second  traverse  mechanism  so  made  that  positioning  of  the  probe  rela¬ 
tive  to  the  ground  board  could  be  made  with  accuracy  (Figure  5). 

The  ground-board  static  pressure  lines  were  connected  to  an  inclined 
manometer  board.  The  boundary  layer  probes,  of  0.  014-inch  outside 
diameter  hypodermic  tubing  with  a  tip  tapered  to  0.  009  inch,  were 
attached  to  a  small  U-tube  manometer  mounted  on  the  reverse  side  of 
the  ground  board.  Pressure  data  from  the  larger  jet  flow  Pitot-static 
probe  were  read  either  from  a  micromanometer  or  from  the  inclined 
manometer. 

Experimental  Results 

Test  data  were  obtained  at  five  ground-board  locations  corresponding 
to  nozzle  height-to-nozzle  diameter  ratios,  H/D  of  4,  2,  1,  0.  5,  and  0.  25. 
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Data  were  obtained  for  two  mass  flows  corresponding  to  jet  nozzle  veloci¬ 
ties  of  121  ft.  /sec.  and  191  ft.  /sec.  At  each  H/D,  velocity  surveys 
were  made  at  the  jet  nozzle  and  at  various  distances  in  the  jet  above  the 
ground  board.  Static  pressure  distributions  were  measured  along  the 
jet  centerline  between  the  ground  board  and  the  jet  nozzle  and  also 
along  the  ground  board.  Measurements  of  velocity  in  the  flow  along  the 
ground  were  obtained  at  various  radial  distances  from  the  stagnation 
point  between  f  =  6  inches  and.  r  =  48  inches. 

Dynamic  pressure  profiles,  where  dynamic  pressure  is  normalized  by 
the  total  pressure  in  the  jet,  are  shown  in  Figure  6  for  the  jet  exhaust¬ 
ing  into  essentially  free  air.  To  obtain  these  data,  the  ground  board  was 
moved  as  far  as  possible  from  the  nozzle  (H/D  =  7.  3).  The  results 
shown  are  for  Um  =121  fps.  There  was  no  measurable  difference 
between  these  profiles  and  corresponding  data  obtained  for  Um  -  191  fps. 
The  figure  shows  that  the  mixing  region  surrounding  the  potential  core 
of  the  jet  grew  in  width  with  increasing  nondimensional  distance  from 
the  nozzle,  ;t/D.  The  potential  core  vanished  completely  at  t/D  **  4.  7. 

Figure  7  shows  the  nondimensional  dynamic  pressure  profiles  in  the  jet 
obtained  with  the  ground  board  set  at  H/D  =  4.  There  was,  again,  no 
difference  between  profiles  obtained  with  64,  =  121  and  191  fps.  Free- 
jet  dynamic  pressure  profiles  corresponding  to  the  same  values  of  t/D 
are  shown  for  comparison.  Note  that  the  two  sets  of  data  are  identical 
for  measurements  made  at  distances  greater  than  about  two  diameters 
■from  the  ground  (z/D>  2).  Evejn  at  z/D  =  1.  08,  the  impinging  jet  corre¬ 
sponds  closely  to  the  free  jet  except  in  the  vicinity  of  the  jet  centerline. 
Deceleration  of  the  impinging  flow  near  the  jet  centerline  at  small 
values  of  z/D  is.  to  be  expected;  the  symmetry  of  the  flow  requires  that 
there  be  a  stagnation  point  at  r/D  =  z/D  =  0. 

Figure  8  shows  the  velocity  distributions  at  the  nozzle  for  H/D  =  4,  2, 

1,  0.  5  and  0.  25;  the  velocities  are  normalized  by  the  jet  reference 
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velocity,  Ua,  .  The  velocity  distribution  for  H/D  =  2.  0  was  essentially 
identical  to  that  for  H/D  =  4.  0.  The  increasing  effect  of  the  ground 
board  on  the  jet  flow  at  the  nozzle  is  quite  evident  in  the  H/D  =1,  0.  5, 
and  O'.  25  velocity  distributions.  At  H/D  =1.0  the  velocity  is  no  longer 
uniform  and  is  reduced  at  the  center  of  the  nozzle  to  0.  96  U <*,  .  The 
velocity  distribution  is  increasingly  nonuniform  for  H/D  =  0.  5  and  0.  25; 
at  H/D  =  0.  25  the  velocity  at  the  center  is  down  to  0.  35  Uao  . 

The  static  pressure  distributions  on  the  ground  are  presented  in  Figures 
9(a)  through  9(f)  for  the  different  values  of  H/D  tested.  Static  pressures 
are  here  made  nondimensional  by  the  static  pressure  difference  between 
the  wind  tunnel  settling  chamber  and  atmospheric  pressure,  which  is 
equivalent  to  the  total  pressure  in  the  jet.  The  data  at  each  H/D  for 
tests  at  U„  =  191  and  121  fps  are  nearly  identical;  hence,  these  pressure 
distributions  are  independent  of  mass  flow.  Although  not  all  points  are 
plotted,  the  data  shown  represent  the  maximum  amount  of  scatter  between 
measurements  made  along  two  perpendicular  lines  through  the  stagnation 
point  and,  also,  between  measurements  made  at  six  evenly  spaced  pres¬ 
sure  taps  on  6-inch  and  12-inch  radii  centered  on  the  stagnation  point. 

The  excellent  symmetry  of  the  flow  in  each  case  is  demonstrated  by 
the  very  small  scatter. 

Figure  10  shows  static  pressure  distributions  along  the  jet  centerline. 

It  can  be  seen  that  the  data,  up  to  H/D  =  2,  lie  along  one  curve,  at  least 
up  to  z/D' s  fairly  close  to  the  nozzle.  The  pressure  distribution  for 
H/D  =  4  is  somewhat  below  the  curve.  Static  pressure  is  essentially 
the  ambient  pressure  for  z/D*Z  1.  5  at  H/D^-2.  Note  that  pressure  data 
were  obtained  a  short  distance  inside  the  nozzle  for  H/D  =  0.  25  and 
H/D  =  0.  50. 

Figure  11  shows  nondimensionalized  velocity  distributions  derived  from 
the  experimental  static  pressure  distributions  measured  on  the  ground 
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board  for  various  values  of  H/D.  These  velocity  distributions  were  com¬ 
puted  by  assuming  that  there  were  no  losses  in  the  flow.  The  velocity 
distributions  for  H/D  =1.0  and  2.  0  are  identical  within  experimental 
error,  while  for  H/D  =  0.  5  and  0.  25  the  distributions  show  the  increas¬ 
ing  effect  of  proximity  of  the  nozzle  to  the  ground. 

Figures  12,  13,  14,  and  15  present  the  boundary  layer  velocity  profiles 
obtained  at  several  radial  stations  between  r/D  =  0.  5  and  r/D  =  1.  33 
for  H/D  =  0.  5,  1,  2,  and  4,  respectively.  No  measurements  were  made 
in  this  range  of  r/D  for  H/D  =  0.  25.  The  velocity  has  been  divided  by 
the  maximum  velocity  in  the  flow  above  the  boundary  layer  at  each  sta¬ 
tion.  The  height  above  the  board  has  been  nondimensionalized  by  the 
nozzle  diameter  and  multiplied  by  the  square  root  of  the  nozzle  Reynolds 
number.  Data  for  both  Ua,  =  191  and  £4e>-  =  121  fps  are  included.  Up  to 
r/D»0.  8,  in  all  cases,  the  boundary  layers  app.ear>~to- 
the  shape  of  the  profiles  and  the  fact  that  data  for  both  mass  flows 
normalize  to  one  curve  with  Reynolds  number.  The  change  in  shape  of 
the  boundary  layer  velocity  profiles  between  r/D  =■  0.  83  and  r/D  =  1.  17 
suggests  that  transition  to  turbulent  flow  has  occurred  in  the  intervening 
region,  although  not  necessarily  at  exactly  the  same  radius  for  each 
value  of  H/D.  Note  the  excellent  repeatability  demonstrated  by  the  data 
of  Figures  12(f)  and  13(f)  and  the  slightly  poorer,  although  acceptable, 
repeatability  of  the  data  in  Figures  14(f)  and  15(f).  The  repeat  data 
were  obtained  after  the  ground  board  had  been  moved  and  then  reposi¬ 
tioned  to  the  H/D  desired. 

The  differences  between  original  and  repeat  data  in  Figures  14(f)  and 
15(f)  can  be  attributed  to  a  change  in  displacement  effect  of  the  boundary 
layer  probe  near  the  ground  plane;  for  both  cases  these  differences  were 
of  the  order  of  0.  001  inch  in  terms  of  distance  from  the  ground  plane. 
Apparently,  at  the  start  of  the  tests,  the  boundary  layer  probe  had  an 
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inner  coating  of  oil  which,  as  the  tests  proceeded,  collected  lint  and 
dust;  as  a  result,  about  halfway  through  the  tests,  the  probe  became 
completely  blocked.  The  probe  was  then  cleaned  out  with  alcohol  and 
a  fine  wire,  and  tests  were  resumed.  The  tests  which  had  been  com¬ 
pleted  prior  to  the  cleaning  were  those  at  H/D  =  4  and  H/D  =  2.  All 
data  for  H/D  =1,  0.  5,  and  0.  25,  and  all  repeat  data,  were  obtained 
after  the  cleaning.  It  is  believed  that  the  lint  and  dust  were  deposited 
asymmetrically  on  the  inside  of  the  probe  in  such  a  manner  that  the 
effective  center  of  the  opening  in  the  probe  was  moved  closer  to  the 
board  than  the  geometric  center.  Consequently,  the  velocity  profiles 
calculated  from  the  data  with  the  geometric  center  of  the  probe  as  a 
reference  would  necessarily  be  displaced  away  from  the  board.  In 
confirmation,  the  theoretical  laminar  boundary  layer  profiles  at  r/D  = 

0.  5  and  0.  667  check  quite  well  with  the  experiments  for  tests  run  after 
the  probe  was  cleaned  and  are  displaced  below  the  experimental  profiles 
by  about  the  same  amount  as  the  difference  between  the  original  and 
repeat  runs  at  H/D  =  4  and  2. 

In  Figures  16(a)  through  16(e),  ujun 7  is  plotted  versus  y/yj£  for  the 
flow  near  the  ground  for  r/D  ^  1.  17  and  -  191  fps  at  all  values  of 
H/D  tested;  yfi  is  the  distance  to  the  point  outside  the  boundary  layer 
at  which  (see  Figure  1).  In  general,  the  data  falls  essentially  on 

a  single  curve  for  1.  33  4  r/D  4  4.  0.  At  H/D  =  4,  the  lower  limit  becomes 
1.  674  r/D.  The  accuracy  of  measurements  made  at  r/D  =  1.  67  is 
doubtful.  To  obtain  the  data  at  r/D  =  1.  67  required  a  modification  of  the 
traversing  mechanism  which  decreased  its  rigidity. 

In  Figures  17(a)  through  17(e),  the  boundary  layer  region  for  the  data  of 
Figures  16  is  plotted  on  a  greatly  expanded  scale.  The  comments  made 
above  apply  here  as  well.  Velocity  profiles  in  this  range  of  r/D  were 
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also  measured  with  =  121  fps.  The  data  obtained  indicated  that  these 
profiles,  when  normalized  as  in  Figures  16  and  17,  are  independent  of 
Uqo  ,  at  least  over  the  range  tested.  A  typical  set  of  data  is  presented 
for  Uco  =  121  fps  and  H/D  =  0.  5  in  Figures  18(a)  and  18(b).  These  figures 
can  be  compared  with  the  corresponding  profiles  for  Uco  -  191  fps. 
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HI  THEORETICAL  ANALYSES 


Basic  Approach 

The  classical  approach  to  the  theoretical  analysis  of  boundary  layers  is 
first  to  compute  the  outer  inviscid  flow  ignoring  the  effect  of  the  boundary 
layer  on  that  flow.  The  boundary  layer  is  then  computed  assuming  that 
the  flow  conditions  on  the  outer  edge  of  the  boundary  layer  are  those  pre¬ 
dicted  at  the  surface  by  the  inviscid  flow  analysis.  This  is  basically 
the  approach  adopted  in  the  theory  reported  in  this  paper.  It  is  assumed 
that  the  effects  of  viscous  mixing  in  the  impinging  jet  flow  do  not  mate¬ 
rially  affect  the  static  pressure  distribution  on  the  ground  plane.  The 
boundary  layer  in  the  outer  ground  flow  region,  where  the  ground-plane 
static  pressure  is  essentially  zero  and  where  viscous  mixing  is  a  basic 
part  of  the  flow  processes,  is  treated  separately. 

Uniform  Inviscid  Impinging  Jet  Theory 

An  inviscid  jet  in  fluid  otherwise  at  rest  is  in  the  class  of  free  stream¬ 
line  flows.  The  essential  difficulty  in  treating  such  flows  analytically 
lies  in  the  fact  that  even  though  there  are  well  defined  boundary  condi¬ 
tions  on  the  free  streamlines  in  such  flows,  the  location  of  the  free 
streamline  is  not  known  in  advance  and  mu^t  be  determined  as  part  of 
the  solution. 

There  is  a  considerable  body  of  literature  concerning  the  theory  of  invis¬ 
cid  and  viscous  jet  flows,  as  noted  in  Reference  5.  Exact  solutions  for 
two-dimensional  inviscid  impinging  jet  flows  using  hodograph  techniques 
have  been  obtained  (References  6  and  7);  but  the  available  theoretical 
treatments  of  inviscid  uniform  three-dimensional  impinging  jets  are 
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limited  to  approximate  or  numerical  methods.  Of  these,  the  approaches 
of  Schach  (Reference  8),  Leclerc  (Reference  9),  and  Shen  (Reference  10) 
are  worthy  of  note  because  the  flow  model  on  which  they  based  their 
analyses  (Figure  19)  is  a  good  approximation  to  the  infinite  impinging 
jet,  and  they  obtain  results  which  are  in  reasonable  qualitative  agree¬ 
ment  with  available  experimental  data. 

Schach  (Reference  8)  used  an  integral  equation  method  developed  by 
Trefftz  and  solved  numerically  for  the  jet  boundary  (BC  in  Figure  19) 
by  successive  approximations.  A  shape  for  the  unknown  boundary  was 
assumed.  On  the  basis  of  the  known  behavior  of  the  velocity  potential 
on  the  assumed  free  boundary  BC  (because  of  the  boundary  condition 
that  the  velocity  is  a  constant)  and  the  boundary  conditions  on  AF  and 
EF,  it  is  possible  to  compute  numerical  values  along  the  assumed 
boundary  BC  of  a  function  which  is  zero  if  the  boundary  is  the  correct 
one  for  the  flow.  If  this  function  is  not  zero,  a  new  jet  boundary  was 
then  assumed,  and  the  process  was  repeated.  The  primary  difficulty 
with  this  technique  was  that  no  rational  basis  was  derived  or  shown  for 
correcting  the  boundary  shape.  In  Schach's  analysis,  the  ratio^  H/D  = 

1.  5  and  l  -  2  in  Figure  19. 

Leclerc' s  approach  (Reference  9)  is  of  particular  interest.  He  used  an 
electrolytic  analog  to  experimentally  determine  the  shape  of  the  free 
streamline  boundary.  The  internal  flow  was  then  computed  numerically 
using  relaxation  techniques.  If  carefully  done,  one  would  expect  the 
experimental  analog  results  to  be  equivalent  to  an  "exact"  solution  for 
the  assumed  model  of  the  impinging  jet. 

Shen  (Reference  10)  obtained  an  approximate  solution  for  the  model  of 
Figure  19,  using  spherical  coordinates,  by  expanding  the  velocity  poten¬ 
tial  f(n  e)  in  series  involving  spherical  harmonics.  A  finite  number 
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of  terms  in  the  series  is  used;  boundary  conditions  are  not  satisfied 
identically  on  the  boundaries,  but  in  an  average  or  weighted  sense.  For 
example,  along  boundary  BC  in  Figure  19,  the  boundary  condition  is 
(d<f>/dS)BC  !-  Ueo  *  where  3  is  tangent  to  curve  BC.  Shen  satisfies  the 
condition 

£  a si  -  sT** <me 

where  is  a  weighting  function.  If  n  terms  in  the  series  for  are 
considered,  there  are  rt  constants  to  be  determined  and,  hence,  /) 
equations  like  the  one  above  are  required.  As  many  equations  are 
obtained  as  needed  by  considering  different  weighting  functions,  . 

On  boundary  BC,  only  the  condition  that  the  tangential  velocity  equals 
t/m  is  used;  for  a  given  assumed  boundary  BC,  the  solution  obtained  is 
checked  by  computing  the  weighted  velocity  component  normal  to  the 
boundary.  If  this  condition  is  satisfied,  within  the  limits  of  the  method, 
the  solution  has  been  obtained.  If  this  condition  is  not  satisfied,  a  new 
boundary  must  be  assumed  and  the  process  repeated.  Shen  indicates  a 
possible  iterative  technique  for  obtaining  the  "exact"  boundary  for  the 
flared  jet  which  involves  the  solution  of  simultaneous  nonlinear  algebraic 
equations,  but  gives  numerical  results  only  for  the  successive  approxi¬ 
mations  type  of  analysis. 

All  of  the  above  authors  used  trial  and  error  techniques  to  determine 
the  free  streamline  boundaries  of  the  impinging  jet  flow,  although  Shen 
does  suggest  a  possible  iterative  technique.  In  addition,  the  effects  on 
the  jet  flow  of  a  ground  plane  in  proximity  to;  a  jet  exit  cannot  be  deter¬ 
mined  from  these  approaches  without  modifications  in  the  assumed  model. 
As  the  jet  exit  approaches  the  ground  plane,  there  is  a  back-pressure 
which  is  reflected  at  the  exit  in  terms  of  reduced  dynamic  pressure  in 
the  center  of  the  jet,  and  the  jet  exit  velocity  profile  is  no  longer  uni¬ 
form  (Figure  8). 
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The  analysis  of  the  inviscid  impinging  jet  performed  during  the  current 
research  was  motivated  by  a  need  for  a  theory  which  would  account  for 
the  effect  of  the  ground  on  the  exit  velocity  distribution  of  a  jet  nozzle. 

A  systematic  iterative  method  to  determine  the  free  jet  boundary  was 
sought.  A  formulation  of  the  inviscid  jet  flow  has  been  derived  in  terms 
of  a  vortex  sheet  representation  for  the  jet  boundaries,  and  a  plausible 
iterative  technique  was  determined  for  deriving  the  free  jet  boundary. 
These  are  presented  in  Appendix  I.  An  IBM  704  program  has  been 
written  and  is  currently  being  checked  out.  At  this  writing,  conclusive 
final  results  have  not  been  obtained;  however  ,  the  results  so  far  obtained 
are  encouraging.  The  method  of  analysis  derived  during  this  research 
may  well  be  applicable  to  a  much  broader  .£llass„  of  inviscid  three- 
dimensional  free -streamline  flows  than  merely  the  uniform  inviscid 
impinging  jet. 


Certain  basic  properties  of  the  flow  of  an  inviscid  impinging  uniform  jet 
can  be  inferred  without  actually  solving  for  the  flow.  Call  0  the  equiva¬ 
lent  uniform  velocity  in  the  tube.  At  some  distance  along  the  ground 
plane  from  the  stagnation  point,  where  the  pressure  has  fallen  to  near 
ambient  pressure,  for  O^lj  ^  Ys(f~),  where  Ys(t)  is  the  height 

of  the  free  stream  surface  above  the  ground  plane.  Reich  (Reference  11) 
showed  experimentally  that  this  is  true  for  r/:D^l.  5.  Conservation  of 
mass  requires  that  ,  _ 

TcpR  O  **  ZrtprYs  Um 


or 


Ys  (r)  « 


;M'D 

ZrU*, 


Hence,  for  a  given  H/D  (and,  hence,  i//Uto  ),  Ys  varies  hyperbolically 
with  t'  »  This  fact  was  utilized  by  Schach  and  Shen  (References  8  and  10) 
in  their  work,  and  Leclerc  (Reference  9)  used  this  result  as  part  of  a 
check  on  the  accuracy  of  the  results  of  his  electrolytic  analog  experiment. 


CD 
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Boundary  Layer  Analysis 

The  boundary  layer  on  the  ground  plane  created  by  an  impinging  jet  can 
be  divided  into  two  distinct  regions  separated  by  a  transition  zone.  In 
the  vicinity  of  the  stagnation  point,  the  boundary  layer  will  be  laminar 
and  its  properties  will  be  determined  by  the  static  pressure  distribution 
on  the  ground.  At  radial  distances  from  the  stagnation  point  sufficiently 
large  to  approximate  a  zero  static  pressure  gradient  at  the  ground,  the 
flow  will  assume  the  form  of  a  radial  wall  jet  which  has  been  analytically 
described  by  Glauert  (Reference  12).  In  the  following  analysis,  the  two 
regions  have  been  considered  separately,  and  an  assumption  about  the 
loss  in  momentum  flux  in  the  transition  zone  has  been  used  to  estimate 
the  initial  radius  of  the  wall -jet  flow  and  its  properties  at  that  radius. 

Laminar  Boundary  Layer 

An  approximate  method  devised  by  Smith  (Reference  13)  was  used  to 
calculate  the  properties  of  the  laminar  boundary  layer  in  the  vicinity 
of  the  stagnation  point.  This  method  was  selected  because  of  its 
accurate  prediction  of  velocity  profiles  in  regions  of  accelerating  flows, 
such  as  exist  in  the  present  application,  and  because  of  its  simplicity. 
The  experimental  pressure  distributions  on  the  ground  were  used  in 
these  calculations  for  each  H/D.  The  potential  velocity  distributions 
associated  with  the  pressure  distributions  have  been  presented  in  Figure 
11.  Boundary  layer  velocity  profiles  predicted  for  various  values  of 
H/D  are  shown  in  Figures  12  through  15. 

The  location  of  transition  from  laminar  to  turbulent  flow  in  the  boundary 
layer  was  investigated  by  consideration  of  boundary  layer  stability 
theory.  Two-dimensional  stability  criteria  were  applicable  to  this 
work  (Reference  4).  These  criteria  were  used  to  determine  the  radius 
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for  neutral  stability  of  the  laminar  boundary  layer,  (r/D)cr^  »  i*1  the 
following  manner.  The  critical  Reynolds  number,  ,  was  obtained 

from  Reference  14  as  a  function  of  Pohlhausen' s  form  parameter,  _A.  . 

A  relation  between  J\-  and  the  included  wedge  angle  for  flow  over  a 
wedge  was  taken  from  Reference  15.  The  equivalent  wedge  angle,  U-rr>l Un 
and  &*/D  had  already  been  obtained  as  functions  of  r/D  in  the  applica¬ 
tion  of  Smith's  method,  so  these  results  were  combined  to  calculate 
curves  of  nozzle  Reynolds  number  versus  (r/D)cr^t  •  Figure 

20  shows  the  final  result  of  these  calculations  for  various  nozzle-to- 
ground  distances. 


In  view  of  the  large  amount  of  turbulence  in  the  mixing  region  of  the 
impinging  jet,  the  stable  region  in  Figure  20  should  be  interpreted  as 
an  indication  of  the  maximum  radial  distance  that  the  boundary  layer 
will  remain  laminar  for  a  given  nozzle  Reynolds  number.  Also,  because 
of  the  steepness  of  the  curves  in  the  region  r/D<l,  it  can  be  inferred 
that  turbulence  feeding  into  the  laminar  boundary  layer  at  values  of  r/D 
only  slightly  less  than  (r/D)cr^t  will  be  highly  damped.  Hence,  the 
curves  should  provide  a  relatively  good  criterion  for  estimating  the 
radius  of  transition  from  laminar  to  turbulent  boundary  layer  flow. 

This  criterion  should  be  especially  accurate  when  (r/D)cr^  K  I. 


Wall  Jet 

At  values  of  r/D£  2,  the  pressure  gradient  along  the  ground  was  found 
to  be  essentially  zero  for  all  values  of  H/D  tested.  In  this  region, 
Glauert's  analysis  (Refer  ence- 12)  of  the  radial  wall  jet  should  apply. 
The  turbulent  flow  solution  is  the  one  of  interest  for  the  present  inves¬ 
tigation.  The  shape  of  the  theoretical  velocity  profiles  was  prescribed 
by  a  parameter,  at  ,  which  in  turn  depended  weakly  on  the  Reynolds 
number,  ^  ,  where  U-m  is  the  local  peak  velocity  in  the  radial 
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wall  flow  and  $t  is  the  distance  between  the  point  at  which  OL**um  and  the 
point  at  which  U^^Um  (see  Figure  1).  A  free  constant,  K ’/  ,  was  used 
to  determine  a  quantitative  relation  between  and  ct  .  The  value 

of  K,  was  found  to  be  Q-.  013  in  Reference  16.  This  value  of  K.  ,  rather 
than  the  value  /C  =  0.  012  suggested  by  Glauert,  was  used  in  the  current 
analysis. 


It  is  desirable  to  obtain  an  estimate  of  ct  from  the  properties  of  the 
jet  before  impingement.  Because  of  the  weak  dependence  of  ct  on  ^  £, 
this  can  be  accomplished  with  considerable  accuracy  even  if  some  rela¬ 
tively  broad  assumptions  are  made.  It  was  assumed  that  Glauert' s 
theory  applies  at  r/D*  1.  5,  that  Ufo/U^ 0»1  at  this  radius,  and  that  the 
loss  in  radial  momentum  flux  due  to  the  wall  shear  stress  is  negligible. 
On  this  basis  ct  =1.16  for  0.  50  H/D  ^  4.  0  and  cL  -  1.17  for  H/D  = 

0.  25  are  predicted.  These  values  compare  favorably  with  the  value 
ct  =  1.  16  which  was  derived  from  the  experimental  measurements  of 
for  0.  25  H/D>^  4.  0.  A  numerical  example  of  this  calculation 

is  presented  in  Appendix  II. 


Once  the  shape  of  the  velocity  profiles  in  the  wall-jet  portion  of  the  flow 
has  been  determined  by  ct ,  the  functional  dependence  of  and  LCm 

on  P*  is  given  by 


c  f—T  y±=c 

U  —  W  /  o!  .  p?  - 


R 


where  the  exponents  a,  and  b  depend  on  the  value  of  cL  ,  and  the  con¬ 
stants  Cf  and  C&  may  be  functions  of  H/Di  The  determination  of  Cf 
and  Q  forms  the  remaining  portion  of  the  analysis. 


(2) 


As  noted  in  Equation  (1),  the  height  of  the  free-stream  surface  of  the 
inviscid  impinging  jet  above  the  ground  plane  is,  for  r/D^  1.  5, 
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wheire 


0* 


■ver'cLir- 


in  the  nozzle. 


The  radial  momentum  flux  My,  is  given  by 

Mr=2  rz/>r  t/s  (r) —,icf>R2‘Uoa  U  (3) 

provided  that  r/D  is  large  enough  so  that  the  pressure  gradient  at  the 
wall  is  approximately  zero.  This  relation  must  necessarily  be  applicable 
to  the  direct  impact  of  two  equal  viscous  jets,  where  there  is  no  mechanism 
for  a  loss  in  momentum  flux,  even  though  the  mass  flow  may  be  greatly 
increased  between  the  nozzle  and  the  radial  station  in  question.  The  jet 
impinging  on  the  ground,  on  the  other  hand,  will  suffer  a  loss  in  momentum 
flux  due  to  the  presence  of  the  wall  shear  stress,  Vq  .  This  loss  can 
be  calculated  for  the  laminar  portion  of  the  boundary  layer.  With  an 
assumption  about  the  additional  losses  in  the  transition  region,  a  momen¬ 
tum  flux  relation  can  be  derived  to  obtain  one  equation  for  the  propor¬ 
tionality  constants  in  Equation  (2).  This  relation  will  have  the  form  of 
Equation  (3)  with  a  term  added  to  account  for  losses  due  to  the  wall 
shear  stress.  It  can  be  written 

Mr  t  /CU  .  j  (4) 

vt&g/ruS  I  V / 

where  Mr  is  the  radial  momentum  flux  at  radius  and  AMr>  is  the  loss 
in  radial  momentum  flux  up  to  V  . 

An  additional  relation  is  required  to  obtain  the  two  unknowns  in  Equation 
(2).  This  relation  was  obtained  empirically.  It  was  observed  that  experi¬ 
mentally  the  quantity  is  independent  of  the  mass  flow  entrain- 

f  \  UcO 

ment  process  occurring  in  the  jet  before  impingement.  This  result 
suggested  that  the  radial  mass  flow  could  be  combined  with  &IU<x>  to 
obtain  a  universal  curve  of  radial  mass  flow,  Or,  versus  7*  .  Now  the 
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radial  momentum  flux  Mp~U-n \  and  the  mass  flow  • 

If  AMp  is  either  negligible  or  roughly  constant  for  all  H/D's  at  the  start 
of  the  wall-jet  flow,  then  -yS  K/*  jr  and  -  ^ 

(/go  t/« 


/  —  /  ,  \/  ^aa ' '  'J*B  U°°R 

tional  to  (U/l/^  where  /  is  a  constant.  Hence,  the  quantity 


will  be  propor- 


<2* 


* ■/>**!/« 


m 


would  be  independent  of  H/D.  The  collapse  of  the  experimental  data 
with  X  =  0.  65  is  illustrated  in  Figure  21.  This  normalization  appears 
to  be  quite  satisfactory  except,  possibly,  for  the  data  obtained  at  H/D  = 
1.  0.  At  values  of  r/R&4,  the  pressure  gradient  is  not  zero,  so  that 
the  wall-jet  solution  does  not  apply  in  this  region. 


The  variation  of  Qp  with  T>  is  prescribed  by  the  form  parameter  CL  . 
From  Reference  12  the  exponents  cl.  and  6  in  Equation  (2)  are  -1.  075 
and  1.  008,  respectively,  for  ct  =  1.16.  Hence,  Qp ~ 

The  effect  of  varying  nozzle  size  and  £4>  on  the  exponent  should  be  neg¬ 
ligible  over  quite  a  large  range  of  these  parameters  because  of  the 
insensitivity  of  cC  to  nozzle  Reynolds  number.  The  relation 


/CU  0  ,  6  fr_ ) 

njoU'+Rk('u)  '*[/?/ 


(5) 


seems  to  give  the  best  fit  to  the  experimental  data. 


Equations  (4)  and  (5)  provide  a  means  of  finding  the  unknown  constants 
Cf  and  C&  .  The  loss  in  momentum  flux,  ,  was  estimated  by 

assuming  that  the  boundary  layer  is  completely  laminar  up  to  the  transi¬ 
tion  point  predicted  in  Figure  20, and  that  the  loss  in  momentum  flux 
varies  with  T*  in  accordance  with  Glauert's  solution  for  the  radial  wall 
jet  beyond  this  radius. 
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After  computing  that 


£7—=  1*6 


Mr, 


K/>J?Uco  KfR'i/J  \  u.m> 

in  the  wall  jet  for  CL  =  1.  16,  one  can  obtain 


U, 


m 


CL 


3.18 


AM, 


C  Tt/iR^U  J 


(±_fss(  rt  j  o.ii« /r  paw 

\W  (~rJ  W 


where  is  the  loss  in  momentum  flux  in  the  laminar  boundary 

layer  up  to  the  predicted  transition  radius,  .  Both  of  these  quantities 
are  available  from  the  laminar  boundary  layer  calculations.  Hence, 


C,=3J8 


1 


L  itjOR'U^O J 


Hit) 


Finally,  since*  =  2.22  C,  C,  fe) 


Q  \  OSS 

CL 

T\0.W 


0.14*. 


(6) 


we  have 


for  CL  =  1.  16  , 


„  _  0.Z07  ( U  \°,6S 
4  \UJ 

These  relations  for  Cj  and  Cz  were  used  in  conjunction  with  Equation  (2) 
to  predict  the  variation  of  and  jR  with  r/R  in  the  wall- 

jet  region.  The  results  are  presented  in  Figures  22  and  23. 


(7) 


The  numerical  constants  in  these  equations  and  Equation  8  on  the  follow¬ 
ing  page  were  obtained  by  integration  of  the  theoretical  velocity  profiles. 
That  is, 

i  .4  -  Inf  um 

where 

Mrm-&K/?rjT-  u*-d(/=  Z rep (uj\  rj  (~j^ 


--poo 

hpfrU',,  Ueo  R  R 

and 


*=1.16 


npR*-Uee>*' 

CL  =  1.16 


1524, 


7.524 
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Transition  Zone 


In  the  zone  between  the  laminar  boundary  layer  region  and  the  wall-jet 
region,  the  ground  flow  was  assumed  to  follow  Glauert's  velocity  pro¬ 
file  with  the  condition  that  0  =  constant  =  value  of  at 

predicted  transition.  This  assumption  and  the  one  pertaining  to  the  loss 
in  momentum  flux  used  in  the  wall-jet  analysis  allowed  an  estimation 
of  y^  jR  in  the  transition  zone.  With  fit  =  1.  16,  the  result  was 

0  }fr 

IZ 


U  l r 73  — — — —  0  (i- 

'  1-52+  (Jec  (  rc/>Ra’UUoal{rtJ  l  <*mt 


-/ 


(8) 


where  LL m is  the  peak  velocity  at  predicted  transition  (0.  98  4  K-  0.  99 
for  all  values  of  H/D  considered  in  the  present  work).  The  ground  flow 
predicted  in  this  manner  cannot  apply  immediately  downstream  of  the 
transition  radius  because  a  finite  distance  is  required  to  establish  some 
sort  of  equilibrium  between  the  jet  mixing  zone  and  the  ground  boundary 
layer.  It  is  plausible,  however,  that  this  distance  may  be  considerably 
less  than  that  required  to  establish  a  pure  wall-jet  flow.  Hence,  the 
above  relation  is  proposed  to  bridge  a  portion  of  the  gap  between  the 
wall-jet  flow  and  the  laminar  boundary  layer. 
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IV  COMPARISON  OF  THEORY  AND  EXPERIMENT, 
AND  DISCUSSION  OF  RESULTS 


Free-Jet  Flow 


Points  of  constant  dynamic  pressure  in  the  mixing  region  of  the  jet 
exhausting  into  free  air  are  shown  in  Figure  24.  These  points  were 
cross-plotted  from  the  dynamic  pressure  profiles  of  Figure  6.  The 
data  appear  to  be  well  represented  by  straight  lines  radiating  from  a 
virtual  origin,  t0  j D  »  -0.  50  an  d  r0/D  »  1.  04,  except  for  stations  very 
close,  to  the  nozzle  exit.  In  the  range  of  tj0  where  this  radial  straight 
line  representation  is  applicable,  the  dynamic  pressure  distributions 
will  be  similar,  that  is,  will  be  a  unique  function  of 

(r-r0)  /tt-to)  .  As  shown  in  Figure  25,  the  dynamic  pressure  dis¬ 
tributions  in  the  mixing  region  are  at  least  approximately  similar  when 
plotted  versus  (r~T$  /  (t~t^  ,  provided  that  the  distance  from  the  nozzle 
exit  is  greater  than  t/D^I  .  The  function 

provides  an  excellent  empirical  fit  to  the  data.  The  velocity  distribu¬ 
tion  is  given  simply  by  the  squre  root  of  this  function: 


where 


_zr 

<7* 


-5~6 


=  -  e 


t-tc 


+  0.097 


$ 


J~  —  -  0.50 


for 


JL>/ 

D?/ 


Kuethe  (Reference  17)  has  derived  theoretically  the  velocity  distribu¬ 
tion  in  the  mixing  region  surrounding  the  core  of  an  axially  symmetric 
circular  jet.  The  theoretical  constant  velocity  curves  are  compared 
with  the  current  data  in  Figure  26.  The  agreement  between  theory  and 
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experiment  is  satisfactory  except,  possibly,  near  the  outer  edge  of  the 
mixing  region  where  the  velocity  is  almost  zero.  Velocity  distribu¬ 
tions  derived  from  the  theory  should  be  more  accurate  than  the  above 
empirical  formula  in  the  region  0  4^(7  ,  whereas, the  empirical 
expression  provides  a  simple  and  probably  quite  accurate  method  of  finding 


velocity  distributions  when  t. 


Impinging  Jet  Flow 


Although  the  theoretical  analysis  of  the  inviscid  impinging  jet  has,  as 
yet,  not  been  completed,  some  calculations  have  been  made.  The 
results  of  calculations  of  ground-plane  and  jet-centerline  pressure 
distributions  based  on  the  vortex  sheet  model  of  the  present  theory 
and  using  the  curve  for  the  free  jet  boundary  computed  by  Schach  are 
presented  in  Figures  27  and  28.  The  corresponding  pressure  distribu¬ 
tions  computed  by  Schach,  Leclerc,  and  Shen  (References  8,  9,  and  10) 
and  the  experimental  pressure  distributions  for  H/D  =  1  and  2  are 
shown  for  comparison. 


At  distances  greater  than  some  fixed  height,  say  Hj  ,  from  the  plane, 
the  infinite  inviscid  jet  impinging  normally  on  a  ground  plane  has 
essentially  a  constant  diameter  and  uniform  velocity  distribution  with 
the  velocity  equal  to  Ug,  ,  the  free  streamline  velocity  at  the  jet 
boundary.  Hence,  for  a  uniform  inviscid  jet  issuing  from  a  nozzle 
at  height  } Hj  above  the  ground  plane,  there  would  be  no  difference  in 
the  impinging  flow  from  that  of  an  infinite  inviscid  jet.  Schach  (Ref¬ 
erence  18)  found  this  H/D  to  be  approximately  2  from  experiments  on 
water  jets  in  air,  which  is  consistent  with  the  present  results  (Figure 
10).  For  an  actual  jet  of  air  into  air  at  rest,  viscous  mixing  occurs. 
For  such  a  viscous  jet,  at  H/D's  large  enough  that  the  jet  is  a  fully 
developed  turbulent  jet  prior  to  impingement,  it  is  clear  that  all  of 
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the  impinging  flow  is  influenced  by  the  viscous  mixing  in  the  jet,  and 
that  the  concept  of  an  "infinite  inviscid  jet"  no  longer  is  relevant.  For 
H/D«l.  0,  it  is  reasonable  to  expect  that,  in  the  neighborhood  of  the 
stagnation  point,  the  flow  along  the  ground  is  not  materially  influenced 
by  the  mixing.  Therefore,  there  must  be  some  intermediate  H/D  at 
which  the  effect  of  the  mixing  in  the  jet  prior  to  impingement  is  first 
felt  in  the  flow  along  the  ground  plane,  as  reflected  in  the  ground- 
plane  pressure  distribution  near  the  stagnation  point. 

It  is  inferred  from  the  present  data  and  from  the  results  obtained  by 
Leclerc  (Reference  9)  that  the  effects  of  jet  mixing  on  the  pressure 
distribution  on  the  ground  are  important  only  for  H/D>  2.  Since  the 
velocity  at  the  nozzle  was  uniform  and  equal  to  for  H/D^  2,  the 

effect  of  the  ground  on  the  flow  at  the  nozzle  was  negligible  (Figure  8). 
Hence,  any  differences  in  the  pressure  distributions  measured  at 
values  of  H/D  >  2  can  only  be  due  to  mixing  in  the  jet  prior  to  impinge¬ 
ment.  The  measured  pressure  distributions  on  the  ground  were  iden¬ 
tical  for  H/D  =  1  and  2  (Figure  11), and  agreed  extremely  well  with  the 
results  obtained  by  Leclerc  for  an  inviscid  jet  at  essentially  an  infinite 
value  of  H/D  (Figure  27).  The  pressure  distribution  obtained  at 
H/D  =  4,  however,  differed  substantially  from  those  measured  at 
H/D  =  1  and  2,  even  in  the  neighborhood  of  the  stagnation  point.  There 
is  a  possibility  that  the  effects  of  mixing  and  of  ground  proximity  can¬ 
cel  each  other  for  1^  H/D^.2  and,  hence,  lead  to  agreement  with  the 
inviscid  results  of  Leclerc.  However,  since  the  effect  of  ground  prox¬ 
imity  is  negligible  at  the  nozzle  for  H/DJ^  2,  with  or  without  mixing 
(cf.  previous  paragraph),  one  must  conclude  that  the  effect  of  mixing 
in  the  jet  is  negligible  for  H/DX.2.  It  should  be  noted  that  the  effects 
of  mixing  for  H/D  =  4  were  evident  in  the  vicinity  of  the  stagnation 
point  even  though  measurements  showed  that  there  was  still  a  potential 
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core  in  the  jet  just  prior  to  impingement.  (Potential  core  diameter  was 
approximately  0.  35D  at  one  nozzle  diameter  above  the  ground.  ) 

Because  mixing  in  the  jet  prior  to  impingement  affects  the  ground  pres¬ 
sure  distribution  for  H/D>  2,  it  appears  that  theoretical  solutions  for 
the  inviscid  jet  will  be  useful  only  for  H/D;S  2.  Any  theory  for  H/D 
much  larger  than  2  must,  somehow,  account  for  viscous  mixing  in  the 
jet  if  it  is  to  result  in  accurate  predictions  of  the  ground  plane  pressure 
distribution.  The  starting  point  for  such  a  theory  could  be  taken  at 
two  to  three  jet  diameters  above  the  ground,  where  the  velocity  profile 
was  found  to  be  identical  to  that  of  the  corresponding  free  jet  (cf>  Sec¬ 
tion  II  and  Figure  7). 

Ground  Flow 

Comparisons  between  the  theoretical  laminar  boundary  layer  profiles 
and  experimental  boundary  layer  profiles  are  presented  in  Figures  12 
through  15  for  H/D  =  0.  5,  1,  2,  and  4.  It  is  evident  that  the  experi¬ 
mental  results  are  predicted  quite  satisfactorily  for  r/D<  0.  7,  if  allow¬ 
ance  is  made  for  the  uncertainty  in  specifying  the  correction  for  the 
distance  of  the  effective  center  of  the  probe  from  the  ground.  These 
"zero  errors"  have  been  discussed  in  Section  II.  As  the  radial  distances 
increase  above  r/D  =  0.  7,  the  agreement  between  experiment  and  the 
laminar  theory  deteriorates  progressively  from  the  outer  edge  of  the 
boundary  layer  inward.  For  the  particular  case  presented  in  Figure  12 
(H/D  =  0.  5),  the  calculated  critical  radius  for  transition  of  the  boundary 
layer  was  r/D  =  1.  08  for  =  191  fps  and  r/D  =  1.  17  for  64>  =  121  fps. 

This  result  suggests  that  in  the  subcritical  range  of  r/D,  the  irregulari¬ 
ties  at  the  outer  edge  of  the  boundary  layer  are  due  to  the  diffusion  of 
turbulence  from  the  jet  mixing  region  into  the  laminar  boundary  layer. 
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As  r/D  is  increased,  this  turbulence  is  able  to  penetrate  deeper  into 
the  boundary  layer  by  virtue  of  the  decreased  stability  until,  at  the 
critical  radius,  the  boundary  layer  becomes  completely  turbulent. 

The  diffusion  process  was  more  marked  at  the  larger  H/D' s,  where  the 
jet  mixing  region  before  impingement  was  wider  and  the  rate  of  change 
of  pressure  gradient  with  r/D  more  gradual.  In  fact,  at  H/D  =  4 
(Figure  15),  the  boundary  layer  appeared  to  be  almost  completely 
turbulent  at  r/D  =  1.  0,  whereas  the  calculated  critical  radii  occurred 
at  r/D  =1.2  and  1.  36  for  U<t>  =  191  and  121  fps,  respectively. 


The  experimental  velocity  profiles  of  the  flow  near  the  ground  for 
r/D>  1.  17  are  given  in  nondimensional  form  in  Figures  16(a)  through 
16(e)  along  with  the  velocity  profiles  predicted  by  Glauert  (Reference  12) 
for  values  of  <£  =  1.1  and  1.  2.  If  it  is  recalled  that  the  appropriate 
value  of  a  was  found  to  be  1.  16,  the  agreement  between  theory  and 
experiment  can  be  seen  to  be  generally  very  good  for  1.  33  <  r/D<  4.  0, 
with  the  lower  limit  increasing  to  1.  67  at  H/D  =  4.  In  the  range 
y ly%  >1*2,  the  experimental  data  diverge  from  the  theory.  This  devia¬ 
tion  has  been  noted  by  other  experimenters  (References  16,  19)  and 
is  probably  due  to  errors  incurred  by  angularity  of  the  flow  with  respect 
to  the  probe,  or  to  the  decreasing  accuracy  of  the  measurements  at 
the  very  low  velocities  encountered. 


The  boundary  layer  region  of  Figure  16  is  plotted  in  Figures  17(a) 
through  17(e).  The  theoretical  profile  for  ct  =  1.16  has  been  added  to 
these  figures.  The  theory  again  appears  to  give  an  excellent  predic¬ 
tion  of  the  flow  for  r/D  >  1.  33,  except  in  the  region  very  close  to  the 
wall.  Near  the  wall  the  turbulent  flow  will  undergo  transition  to  what 
is  generally  referred  to  as  the  laminar  sublayer.  The  outer  edge  of 
this  transition  region  can  be  estimated  from  the  criterion  (y^ZQjp)J^  >70 
for  purely  turbulent  flow  (Reference  14).  With  the  aid  of  Blasius' 
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formula  for  turbulent  pipe  flow,  this  becomes 


y  >  w 


■frM'/w 


%  io.ozzs 

for  the  region  of  completely  turbulent  flow  in  Figures  17(a)  through 
17(e).  This  expression  is  approximately  equal  to  0.  015  for  all  values  of 
H/D  tested  in  the  current  work.  Hence,  the  experimental  data  can  be 
expected  to  diverge  from  that  predicted  by  the  wall -jet  theory  when 

0.  015.  Accordingly,  the  theoretical  profile  for  cL  -  1.  16  has 
been  modified  by  matching  the  inner  portion  of  the  universal  velocity 
distribution  for  smooth  pipes  (Figure  20.4  of  Reference  14)  to  the 
theoretical  curve  at  y ) =  0.  015.  The  resulting  curve  now  agrees 
with  the  experimental  data  even  in  the  region  where  y j y-%  is  small 
enough  to  make  the  flow  primarily  dependent  on  viscous  stresses  rather 
than  Reynolds  stresses. 


The  shape  of  the  nondimensional  velocity  profiles  has  been  shown  in 
the  previous  paragraphs  to  be  adequately  predicted  by  theory.  The  ade 
quacy  of  the  analysis  is  also  dependent  on  its  accuracy  in  predicting 
Ufnf  U<b  and  yjg  /R  in  the  wall -jet  region  and  on  whether  Um  is 
equal  to  LAj  in  the  laminar  region.  The  theoretical  variation  of 
Ufn/ Uao  with  r/D  for  Uoa  =  191  fps  is  compared  to  its  experimental 
equivalent  in  Figures  22(a)  through  22(e).  In  the  laminar  region  the 
assumption  (Am**  £2/  appears  to  be  valid  up  to  the  predicted  transition 
radius  for  H/D  =  0.  5  and  1.0.  At  H/D  =  2  and  4,  the  experimental 
value  of  lAn j  begins  to  fall  below  U,j  at  radii  less  than  that  for  pre¬ 
dicted  transition.  This  behavior  is  consistent  with  the  appearance  of 
an  apparently  forced  transition  of  the  laminar  boundary  layer  at  radii 
less  than  that  for  neutral  stability.  In  the  wall-jet  region,  LCfn  j  Uao 
was  evaluated  from  Equations  (2)  and  (6).  The  calculated  curves  of 
um  / Uqo  are,  in  general,  about  10  percent  higher  than  the  experimental 
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data.  Judging  from  these  results,  the  assumption  that  the  momentum 
flux  variation  with  V  has  the  same  form  throughout  both  the  transition 
region  and  the  wall-jet  region  may  possibly  be  inadequate.  However, 
the  errors  could  also  be  due,  in  part,  to  the  occurrence  of  the  above- 
mentioned  forced  transition  or  to  small  differences  between  theory  and 
experiment  in  the  values  of  the  exponents  Cl  and  6  in  Equation  (2). 
Recalculation  of  Cj  with  q,  =-l.  143  and  b  -  1.  028  to  conform  with 
the  experimental  results  provided  estimates  of  Cf  which  were,  in 
general,  considerably  more  accurate  over  the  H/D  range  tested.  The 
curves  obtained  with  this  modified  analysis  are  also  shown  in  Figures 
22.  The  improved  accuracy  suggests  that  the  discrepancy  between 
the  analysis  and  experiment  is  due  primarily  to  small  differences 
between  predicted  and  actual  values  of  a,  and  b  . 


The  variation  of  tR  with  r/D  in  the  wall -jet  region  is  shown  in 
Figures  23(a)  through  23(e).  By  virtue  of  the  inverse  relation  between 
Cf  and  Cz  in  Equation  (7),  these  predictions  are  too  low  by  the  same 
amount  that  u.^1  U ^  is  too  high.  The  improvement  in  the  analysis 
obtained  by  using  the  experimental  values  of  Cl  and  b  is  reflected 
here  as  well  as  in  the  calculation  of  / Uoo • 


In  deriving  Equation  (8)  for  the  value  of  y^/R  in  the  transition 
zone,  it  was  assumed  that  Glauert's  theoretical  velocity  profiles  applied, 
even  though  the  flow  in  this  region  does  not  actually  conform  to  a  true 
radial  wall  jet.  Some  insight  as  to  the  accuracy  of  this  assumption 
can  be  obtained  by  inspection  of  the  velocity  profiles  in  Figures  16  and 
17.  At  all  values  of  H/D  except  4,  the  profiles  are  accurately  predicted 
for  r/D>  1.  33.  AtH/D  =  4  the  boundary  layer  region  of  the  velocity 
profiles  is  predicted  satisfactorily  for  r/D^  1.  33  and  the  outer  flow 
region  is  given  quite  well  for  r/D>  1.  67.  The  value  of  OLm  /  Uoo 
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assumed  for  the  transition  region  is  shown  for  each  H/D  in  Figures 
22(a)  through  22(e),  and  the  resulting  values  of  are  shown  in 

Figures  23(a)  through  23(e).  The  over-all  accuracy  of  the  predictions 
in  this  region  is  comparable  to  that  for  the  wall-jet  flow  provided  that 
r/D>  1.  33.  The  juncture  between  the  theoretical  transition  and  wall- 
jet  regimes  appears  to  give  a  reliable  estimate  of  the  radius  at  which 
the  experimental  data  approach  the  equivalent  of  a  true  wall  jet.  The 
region  between  r/D  =  1.  33  and  the  transition  radius  cannot  be  analyzed 
by  the  method  presented.  Fortunately,  this  region  appears  to  be  very 
small. 
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APPENDIX  I 


DERIVATION  OF  EQUATIONS  FOR  INV1SCID  IMPINGING  JET  FLOW 


Basic  Formulation 


Consider  a  Cartesian  coordinate  system  with  the  z-axis  directed  along 
the  axis  of  symmetry  of  the  flow,  positive  vertically  up  from  the  ground 
plane,  and  the  x-y  plane  (z  =  0)  coincident  with  the  ground  plane.  Con¬ 
sider  an  element  of  volume  <£z  with  vorticity  vector  £  fa '\y[z)  . 

For  vorticity  distributed  throughout  a  given  volume  /  ,  one  can  write 
the  vector  velocity,  ^  ,  at  any  point  (x,  y,  z)  in  the  fluid  as  (Reference 

6) 


where 


dvfx,  y\z.') 


d-1) 


d-  (%-x/)  t  +  (y-y')j  +  (z  ~z')  & 

and  A  are  unit  vectors  in  the  x,  y,  and  z  directions,  respectively. 

In  the  case  of  axially  symmetric  flow. 

The  transformation  to  cylindrical  coordinates  (/]  Zt  B) ,  where 


and 


gives 


and 


x  -  r  cos  6,  y  -  r  sin  Q 

'C0  C«S  0 

£>xd  -  Co\_-*r  (z~z>)  cos  &  ~J9  (z~z')  cos  d  +  (rcos 

I d\3=  [(zrzT  +rz  +  (r')z-  2  rr' cos  0']* 
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where  JB  and  are  unit  vectors  in  the  r  >  &  and  z  directions, 
respectively.  Because  of  axial  symmetry,  only  the  0-0  plane  need 
be  considered. 


Substituting  in  Equation  (I- 1 )  and  equating  vector  components,  we  have 

H_  f  C  (z-z')  cos  Q'd  v  (r,'z'  O') 

4nJv  l(z-z.')A+  r*+  (r')z  -  Zrr' cos&']y* 

1  r  t0  (r,'z')(r  cos  Q'-r')  dr(r'z\  O')  t1"2) 

4tcj^  [{z-z')*'+  r*-+(rj1-. 2  rr'cos  O'] ^ 

O 

where  LL.  ,  'O'  are  the  velocity  components  parallel  to  the  r-axis  and 
the  z~axis,  respectively,  and  t/'q  is  the  velocity  component  perpendicu¬ 
lar  to  the  meridian  plane.  The  integration  is  taken  over  the  volume 
V  which  includes  all  the  distributed  vorticity  of  intensity  £  faz)* 


Consider  now  the  vorticity  distributed  in  a  thin  layer  on  the  edge  of  the 
jet.  Consider  a  given  element  of  volume,  d'C  ,  of  this  layer  (Figure 
29)  whose  projection  on  the  meridian  plane  is  eds  ,  where  £  is  the 
infinitesimal  thickness  of  the  layer  and  ds  is  the  length  of  the  element 
along  the  curve  r  —  P0  (z)  which  defines  the  edge  of  the  jet.  From 
Stoke's  Theorem  (Reference  6,  page  46)  for  a  circuit  C  in  the  meridian 
plane  enclosing  the  projection  of  oLV  ,  we  have 


f'li  ■  f  dA  =  J  <i  •  ds 


where  Fj  is  the  unit  outward  normal  vector  to  the  surface  enclosed  by 
C,  S'  is  the  area  enclosed  by  C,  and  is  the  velocity  vector  around 
C.  The  integral  on  the  right-hand  side  is  taken  in  the  clockwise  direc¬ 
tion.  Then 


/  n  't,cLA  =  C0  eds 
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hence, 


Z°£  "i,~% 

where  ^  and  ^  are  as  shown  in  Figure  29.  If  we  let  €  -*  0»  C0— *  00 
such  that 

%06  -  % •  a  constant 
then  i. 

W,-«. 

A1SO’  C0€r0dsdO' 

-  Zr0(z')-/l  +  [^2 % dz'd.  6 


Thus,  Equation  (1-2)  can  be  written 


u.  (r,z)  =  —~i 
4n 

JS 

v-(r,z)=~ 


Z(r')r0(z‘)  -j] + 

’  .  i 

« _ 1 

iz-z')  cosG'cLz'dO' 

\lz-z')z  i-rxf 

r0  (z ')}z-2rr0  (z‘)  cos  O' 

r 

(1-3) 

iz'd& 

t,(z')r0(z')-/l  + 

‘dr0(z’)' 
■  dz 

*  [r  cos  Q'-r0  (z')]  < 

[(z -z'Y  +  rx  -f-[r0(z')}z- 2 rr0  (z')cos  O'] 

% 

where  the  integration  is  carried  out  over  the  entire  vortex  sheet  *5  . 


The  Boundary  Conditions  for  Infinite  Impinging  Jet 

The  boundary  conditions  on  the  inviscid  impinging  jet  flow  are:  along 
the  surface  z  =  0,  0;  along  the  jet  boundary  which  is  a  free  stream¬ 

line,  the  velocity  is  equal  to  a  constant  and,  therefore,  £  is  also  con¬ 
stant,  as  just  outside  the  jet  the  velocity  is  zero. 

To  satisfy  the  boundary  condition  that  2/^=0  for  z  =  0,  the  plane  z  =  0 
is  made  a  plane  of  symmetry  for  the  impinging  jet  and  its  image.  Hence, 
the  integration  over  S  ,  Equation  (1-3),  must  be  carried  out  over  two 


37 


63-11 


vortex  sheets,  that  defining  the  imping  jet  and  that  for  the  image  jet. 


The  Boundary  Conditions  for  Jet  Nozzle  at  Finite  Height  Above  the  Gr'oUhd' 

The  jet  is  assumed  to  issue  from  a  straight  tube  of  constant  diameter 
(Figure  30)  at  finite  height,.  H  ,  above  the  ground.  A  free-stream  sur¬ 
face  issues  from  the  edge  of  the  jet  exit,  B  in  Figure  30,  such  that  the 
free  streamline  velocity  vector  at  B  is  parallel  to  the  side  of  the  tube. 

The  flow  inside  the  tube  along  the  side  of  the  tube  must  be  tangential 
to  the  side;  but,  of  course,  the  streamline  along  this  side  is  no  longer 
a  free  streamline,  as  a  pressure  difference  can  exist  and,  hence,  the 
magnitude  of  the  edge  velocity  may  vary.  The  boundary  condition  along 
the  edge  of  the  tube  is,  then,  that  the  normal  component  of  velocity, 

U  ,  is  zero.  In  addition,  at  B  in  Figure  30,  V"  is  continuous.  Again, 
along  the  free  streamline  BC,the  velocity  and,  hence,  £  are  constant. 


Final  Form  for  Velocity  Components 


Figure  30  illustrates  the  flow  model  as  thus  far  developed.  The  straight 

tube  of  constant  diameter,  D,  centered  on  the  z-axis  extends  from  x  =  H 

to  z  =  o o  .  For  0<z<  H ,  the  curve  r—T0(z)  defines  the  edge  of  the 

jet.  To  insure  that  V'  is  tangential  at  B  in  Figure  30,  we  require  that 

-^-2  =  0  at  z  =  H;  that  is,  there  is  no  discontinuity  in  slope  between  the 

oiz 

side  of  the  nozzle  and  the  free-stream  surface.  The  sides  of  the  tube 
and  the  jet  are  represented  by  a  vortex  sheet  of  vorticity  Z(r0,z)  per 
unit  area  whose  r-  and  z-components  are  zero.  The  image  of  this 
vortex  system  for  0>z}-co  assures  that  the  boundary  condition  for  no 
normal  flow  across  OD  is  satisfied. 
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The  integrations  with  respect  to  d'  can  be  done  in  closed  form  in 
each  of  the  integrals  of  Equation  (1-3).  If  we  let  Z  m  z.  R  ,  z' ■*"  Z-' R  , 
rmrR,  r0  **  R  and  £-f  %(H)>  there  results 


uMmSiF 


dLz! 


i 


00  %(z-')(z+z') 

(H  jli+zT+fr+l)* 


(z  f2')V  r,z y  / 


£Y4J  -K(AZ) 


cLz‘ 


rh^Ju[^ 


(l  ^z-z'JVlr  f  r0  ^z'/J*  j 


[z-2'f /  ZV[r^?£ 

[2-Z']V[f-r0SV]‘ 


E(A3)-K(A3))dt 


[z^Z'JV  z-Vfr.S'F 


(it 


ir(Kl)= 


1(H) 

Ztt. 


IM 


Wtt 

R 


1/(2 -Z')z  +■  (F  +  U* 

Ifr) 


(z-z')*’+  r*-f 

ITFFTw  *<»-**-> 


dLz' 


(Z  +  z')z  +  rz-1  , 

)h  t/(z  +  z')*  +  (rtf)*  \(z+  Z')1  +  (f~/)z  ^  ^  z 


<Lz 
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where  K (A)  is  the  complete  elliptic  integral  of  the  first  kind,  and 
E(A)  is  the  complete  elliptic  integral  of  the  second  kind,  with  modulus 
■A  ;  also,  ,» _  4rT  /A 

-  (z  -rY  +  lrt  f)*  ’  W  = 


W  =  JZ 


4rf0  (z>) _  //  _  _ 


lz-r)1  -f  [ftfa(2')T 


(*  +r)z  t  JfTi. r 
4rrJif)  , 


(E-2')v  +  [f>  f0  (2')] 


These  integrals  are  valid  everywhere  except  on  the  boundaries,  at  (/,  z) 
for  2  >  %  and  fr0  (z),  £  for  z.  4.  g  ,  where  the  integrands  are  singular. 
However,  even  on  the  boundaries,  it  can  be  shown  that  the  Cauchy  prin¬ 
cipal  values  of  these  integrals  exist. 


Equations  (1-4)  enable  the  velocity  at  any  point  in  the  flow  to  be  com¬ 
puted,  once  the  form  of  (z)  ,  the  curve  defining  the  boundary  of  the 
jet,  and  £  (z)  ,  the  vorticity  area  density  of  the  vortex  sheet  on  the 

boundary  of  the  jet  and  the  tube,  are  known. 


Equations  (1-4)  with  the  associated  boundary  conditions 


tint 

*->1 

Um 

f-+rQ 


u  (r>  z)  =  0 ,  r<1 

_  Uf0(z)  / 

r(ftz.)  -  dz  '  *.<"//<> 


f  <r0 


give  rise  to  simultaneous  nonlinear  integral  equations  for  £/%)  and 
tg  (z)  ,  which,  as  they  stand,  do  not  admit  of  solution  in  closed  form. 
However,  by  means  of  a  high-speed  large -capacity  digital  computer,  a 
solution  may  be  possible  if  a  convergent  iterative  process  can  be  derived 
which  will  permit  the  determination  of  the  jet  boundary  (z)  •  What 
appears  to  be  a  rational  process  for  doing  this  has  been  derived  and 
programmed  on  the  IBM  704  computer.  However,  the  results  obtained 
as  yet  are  not  conclusive  because  of  difficulties  with  the  program. 
Therefore,  only  a  brief  outline  of  the  computational  process  is  presented 
here. 
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The  general  computational  procedure  adopted  is  as  follows: 

1.  The  initial  shape  of  the  free  streamline  boundary  is  assumed. 

2.  The  boundary  condition  that  the  flow  at  the  wall  of  the  tube 
is  parallel  to  the  tube  is  used  to  determine  the  varying  area 
vorticity  density  of  the  vortex  sheet  representing  the  tube 
and  the  constant  area  vorticity  density  of  the  free-stream 
surface  vortex  sheet. 

3.  The  flow  normal  to  the  free-stream  surface  is  then  evaluated; 
if  there  is  a  velocity  normal  to  the  boundary  assumed  for 

the  free-stream  surface,  that  assumed  boundary  is  incorrect 
and  must  be  adjusted. 

4.  The  iteration  proceeds  for  the  adjusted  boundary  (by  going 
back  to  Step  1). 

The  process  of  determining  the  revised  jet  boundary  in  Step  3  is  based 
on  the  boundary  condition  along  the  jet  boundary  free  streamlines,  that_ 
the  velocity  normal  to  the  boundary  is  zero.  For  example,  if  the  com¬ 
puted  velocity  normal  to  the  jet  boundary  is  directed  outward,  this 
intuitively  would  suggest  that  the  assumed  boundary  should  be  displaced 
outward.  Let  the  curve  for  the  assumed  free  jet  boundary  BC  (Figure 
30)  for  the  i**1  repetition  of  Steps  1  through  3  be 

(rJe  *  [r0&)i 

Once  the  distribution  of  the  vorticity  in  the  vortex  sheet  is  known  from 
Step  2,  above,  it  is  possible  to  compute  the  velocity  vector  anywhere 
in  the  flow,  based  on  (K  k  ,  the  assumed  free-stream  surface  boundary 
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using  Equations  (1-4).  The  velocity  components  u(fot  2 )  ,  V"  (for  Z) 

at  the  boundary  are  computed;  then  \u  (f0  >z) / v-  (re,  z)J  is  the  slope 

tVi 

of  the  velocity  vector  along  the  itn  assumed  boundary.  If  this  slope 
is  equal  to  JcLz.  all  along  the  assumed  boundary,  (f0)i  is  the 

correct  curve,  as  all  boundary  conditions  are  then  satisfied.  If  not,, 
two  possible  approaches  for  deriving  a  revised  boundary  are  as  follows: 


1.  As  r0 


« _  i  .r  ,  o<t<£ 


a  new  boundary  &0)i+f  can  be  obtained  from 

u(re.  ,  2)  _  cL  [r0(z')]t- 
v(rai  ,z)  d.  z 


di' 


where  K \  is  a  constant.  The  slope  of  [f0/zj J.^  is  theii 

d[r„(z)\i¥/  _  (j -  y)  ^  ^ +  %  jj  fc;  ,  z) 

dii  '  di.  1  tz) 

2.  More  simply,  we  could  let 


(1-5) 


d-6) 


where  K ,  is  a  constant. 


Either  approach  (1  or  2),  hopefully,  provides  a  rational  basis  for  an 
iterative  technique.  It  should  be  noted  that  no  theoretical  rigor  is 
claimed  for  this  procedure;  it  is  primarily  a  systemization  of  the  method 
which  a  computer  might  use  to  adjust  boundaries  in  a  successive  approxi¬ 
mation  technique.  It  is  obviously  intended  for  use  in  an  electronic 
digital  computer  program.  Convergence  can,  thus  far,  be  determined 
only  by  a  trial  and  error  process.  Several  iterations  would  be  made 
with  a  selected  value  of  either  of  the  constants  K,  or  Kt  depending  on 
which  method  is  used.  If  the  successive  boundaries  appear  to  be 
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diverging,  then  successively  smaller  values  of  Kt  or  would  be 
used  until  convergence  is  obtained; 

Method  1,  above,  (Equation  1-5)  for  deriving  successive  jet  boundaries 
was  incorporated  in  the  IBM  704  program  with  the  constant  =1.0, 
In  this  program,  the  distribution  of  vorticity  density  £  is  determined, 
for  a  given  assumed  free  jet  boundary,  by  a  collocation  procedure.  At 
specified  points  on  the  tube  (boundary  AB  in  Figure  30),  the  velocity 
normal  to  the  tube  is  computed  in  terms  of  ?  and  set  equal  to  zero, 
There  results  a  set  of  linear  equations  in  terms  of  £  (z) .  As  a  check 
on  the  integration  accuracy,  the  velocity  at  the  free  jet  boundary  is 
computed,  once  convergence  is  obtained,  to  determine  whether  it  is, 
in  fact,  equal  to  Uoo  . 

As  noted  in  the  body  of  this  report,  calculations  of  velocities  along  the 
jet  centerline  and  along  the  ground  plane  have  been  made  using  as  an 
assumed  curve  for  the  jet  boundary,  one  derived  by  Schach  (Reference 
8)  by  a  successive  approximations  technique.  As  this  boundary  corre¬ 
sponds  essentially  to  that  for  the  infinite  impinging  jet,  the  correspond¬ 
ing  vorticity  area  density  on  the  vortex  sheet  should  be  constant,  and 
this  was  assumed.  The  ground  plane  static  pressures  and  jet  center- 
line  static  pressures  resulting  from  this  calculation  are  shown  in 
Figures  27  and  28  and  are  in  relatively  good  agreement  with  experi¬ 
mental  pressures  for  H/D  =  2.  The  agreement  is  better,  generally, 
than  Schach's  computed  pressures.  The  difference  between  the  present 
results  and  those  of  Schach,  based  on  the  same  boundary,  are  probably 
a  result  of  the  more  accurate  numerical  integration  which  can  be 
obtained  on  a  large  digital  computer. 

Preliminary  calculations  obtained  during  partial  runs  of  the  complete 
program  are  encouraging  and  offer  hope  that,  in  fact,  the  process  will 
converge.  It  is  expected  that  eventually  calculations  of  the  inviscid 
impinging  jet  flow  will  be  obtained  for  a  range  of  values  of  H/D. 
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APPENDIX  II 


CALCULATION  OF  GLAUERT'S  FORM  PARAMETER,  CL 


The  value  of  the  shape  parameter  dL  as  a  function  of  the  Reynolds 
number  LLnfS^j-^  is  shown  in  Figure  31.  This  curve  was  calculated 
from  the  data  presented  in  Table  1  of  Reference  12. for  K*  =  0.  013. 

It  is  evident  that  cC  depends  only  weakly  on  itfn&ij'P  .  Thus,  only  a 
rough  estimate  of  the  Reynolds  number  is  required  if  cL  is  to  be 
determined  with  satisfactory  accuracy;  an  error  of  -0.  02  in  at  will 
lead  to  negligible  differences  in  the  predicted  characteristics  of  the 
flow  at  the  ground. 


It  has  been  generally  observed  that  the  flow  on  the  ground  under  an 
impinging  jet  reaches  a  maximum  velocity  at  IX  r/D&2  and  that  this 
maximum  velocity  is  approximately  equal  to  that  corresponding  to  the 
total  pressure  in  the  jet  before  impingement.  With  these  observations  as 
a  guide,  it  was  assumed  that  U^jUoa  -  1  at  r/D  =  1.  5,  that  Glauert's 
theory  applies  at  this  radius,  and  that  the  loss  in  momentum  flux  due 
to  the  wall  shear  stress  is  negligible  up  to  this  radius. 


The  radial  momentum  flux  inwall-jet  flows  is  given  by 

Mr 


-  faff  £  ■ 


7t/>R*Utxsl 

The  constant  is  equal  to  1.  514  at  oi  -  1.0  and  increases  gradually  to 
1.  547  at  at  -  1.  6.  Since  this  range  of  at  is  the  one  most  likely  to  be 
encountered  in  practice,  an  average  value  of  the  constant,  equal  to 
1.  53,  represents  a  very  good  compromise.  Hence, 

Mr 


7tf>RxUi 


CO 


m  f  r 


44 


63-11 


and  from  Equation  (4)  in  the  body  of  this  report 


u 


64, 

Combining  these  two  equations  and  putting 


if  losses  are  neglected. 


jj0-  -  1  at  r/D  =  1.  5  gives 

C/oo 


“my* 

U^R 


2* 


153 


'</*> 


or 


&t  _  y^/t. 

S  7* 

A  first  approximation,  (at)f 


d)  Its) 

&)  *  m  (£)(Y)(ir)  ■ 


CL 


to  the  required  value  of  cl  is  obtained 
by  choosing  Sf -  1  and  using  Equation  (II - 1 )  to  find  —  jfg. /*/&■ 

and,  hence,  (cL ,)  from  Figure  31.  For  example,  in  the  current  tests, 
at  H/D  =  0.  5  and  Z/«  =  191  fps,  J&l£  equalled  1.  22  x  106  and  O/i/oo 

was  0.84.  This  gave  (/— ^r~)  =  1.  12  x  10^,  which  corresponds  to  d  = 

1.  16.  At  this  value  of  (cL)f  ,  Stj «*  0.  9  and  a  second  approximation 
using  £?r£  =  »•  ’  can  be  performed.  The  result,  (*h  - 

1.  16,  is  the  same  as  that  for  (cC)t  .  At  very  much  lower  Reynolds 
numbers,  the  second  approximation  will  be  slightly  higher  than  the 
first. 


At  very  large  values  of  H/D,  the  jet  just  prior  to  impingement  will  be 
fully  turbulent.  In  this  case,  effective  values  for  nozzle  radius  and 
C/oo  should  be  used.  The  effective  radius  could  be  defined  by  the  radius 
of  a  free  jet,  at  a  distance  from  the  nozzle  equal  to  H,  for  which  the 
velocity  equals  j  the  centerline  velocity.  The  reference  velocity  Uoo 
would  be  taken  as  the  free  jet  centerline  velocity  at  this  same  distance. 
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STATIC  PRESSURE 
DISTRIBUTION 


Figure  l.  NOTATION  FOR  FLOW  IN  IMPINGING  JET 
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Figure  2.  VIEW  OF  EXPERIMENTAL  APPARATUS 
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Figure  3.  VIEW  OF  EXPERIMENTAL  APPARATUS 


Figure  5.  BOUNDARY  LAYER  TRAVERSE  MECHANISM 


Figure  7.  DYNAMIC  PRESSURE  PROFILES  IN  IMPINGING  JET  (H/D  =  4.0) 
COMPARED  WITH  FREE- JET  PROFILES 


Figure  9(d).  GROUND-PLANE  STATIC  PRESSURE  DISTRIBUTION  -  H/D  =  2.0 


■re 
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Figure  12(a).  GROUND-PLANE  BOUNDARY  LAYER  Figure  12(b).  GROUND-PLANE  BOUNDARY  LAYER 
PROFILES  -  H/D  =  0.50,  r/D  =  0.5  PROFILES  -  H/D  =  0.5,  r/D  =  0.667 


Figure  13(a).  GROUND-PLANE  BOUNDARY  LAYER  Figure  13(bX  GROUND-PLANE  BOUNDARY  LAYER 
PROFILES  -  H/D  =  1.0,  r/D  =  0.5  PROFILES  -  H/D  =  1.0,  r/D  =  0.667 


Figure  13(c).  GROUND-PLANE  BOUNDARY  LAYER  Figure  13(d).  GROUND-PLANE  BOUNDARY  LAYER 
PROFILES  -  H/D  =  1.0,  r/D  =  0.83  PROFILES  -  H/D  =  1.0,  r/D  =  1.0 


Figure  14(c).  GROUND-PLANE  BOUNDARY  LAYER  Figure  14(d).  GROUND-PLANE  BOUNDARY  LAYER 
PROFILES  -  H/D  =  2.0,  r/D  =  0.83  PROFILES  -  H/D  =  2.0,  r/D  =  1.0 
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Figure  15(o).  GROUND-PLANE  BOUNDARY  LAYER  Figure  15(b).  GROUND-PLANE  BOUNDARY  LAYER 
PROFILES  -  H/D  =  4.0 ,  r/D  =  0.5  PROFILES  -  H/D  =  4.0,  r/D  =  0.667 


>|<5> 


3 


71 


63-11 


Figure  15{c)  GROUND-PLANE  BOUNDARY  LAYER  Figure  15(d).  GROUND-PLANE  BOUNDARY  LAYER 
PROFILES  -  H/D  =  4.0,  r/D  =  0  .83  PROFILES  -  H/D  =  4.0,  r/D  =  1.0 


Figure  15(e).  GROUND-PLANE  BOUNDARY  LAYER  Figure  15(f).  GROUND-PLANE  BOUNDARY  LAYER 

PROFILES  -  H/D  =  4.0,  r/D  =  1.17  PROFILES  -  H/D  =  4.0,  r/D  =  1.33 


Figure  16(a).  VELOCITY  PROFILES  NEAR  THE  GROUND  FOR  r/D  >1.0 

H/D  =  0.25,  Uw  =  191  FPS 
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igure  16(b).  VELOCITY  PROFILES  HEAR  THE  GROUND  FOR  r/D  >1.0 

H/D  =  0.50,  U*  =  191  FPS 


Figure  16(e).  VELOCITY  PROFILES  HEAR  THE  GROUND  FOR  r/D  >1.0 

H/D  =  1J) ,  U.  =  191  FPS 


sad  i6i  =  °°n  '<rz  =  a/H 


Figure  17(b).  BOUNDARY  LAYER  VELOCITY  PROFILES  FOR  r/D>1.0 

H/D  =  0.50,  U„  =  191  FPS 


H/D  =  1.0  .  U*  =  191  FPS 


Figure  17(d).  BOUNDARY  LAYER  VELOCITY  PROFILES  FOR  r/fo>1.0 

H/D  =  2.0  ,  Uoo  =  191  FPS 


Figure  17(e).  BOUNDARY  LAYER  VELOCITY  PROFILES  FOR  r/D»1.0 


Figu*  18(a).  VELOCITY  PROFILES  HEAR  THE  GROUND  FOR  r/D  >1.0 

H/D  =  030.  U«,  =  121  FPS 


H/D  =  050.  Uo.  =  121  FPS 


UA 


RADIAL  POSITION  OF  THE  POINT  OF  NEUTRAL  STABILITY 
IN  THE  LAMINAR  BOUNDARY  LAYER  VERSUS  NOZZLE 
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Figure  21.  UNIVERSAL  CURVE  OF  RADIAL  MASS  FLOW  vs  r/R 
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Figure  22(o).  VARIATION  OF  EXPERIMENTAL  MAXIMUM  LOCAL  VELOCITY  NEAR  THE 

GROUND  WITH  r/D  COMPARED  WITH  THEORY  FOR  H/D  =  0.25,  U«  =  191  FPS 
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Figure  22(b).  VARIATION  OF  EXPERIMENTAL  MAXIMUM  LOCAL  VELOCITY  NEAR  THE 

GROUND  WITH  r/D  COMPARED  WITH  THEORY  FOR  H/D  =  0.5  ,  U«  =  191  FPS 
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Fiflw*  22(d).  VARIATION  OF  EXPERIMENTAL  MAXIMUM  LOCAL  VELOCITY  NEAR  THE 

GROUND  WITH  r/D  COMPARED  WITH  THEORY  FOR  H/D  =  2.0  ,  Uoo  =  191  FPS 
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Figw*  22(«X  VARIATION  OF  EXPERIMENTAL  MAXIMUM  LOCAL  VELOCITY  NEAR  THE 

GROUND  WITH  r/D  COMPARED  WITH  THEORY  FOR  H/D  =  4.0  ,  Uo*  =  191  FPS 
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Figure  23(a).  VARIATION  OF  EXPERIMENTAL  VALUES  OF  y^/R  WITH  r/t) 
COMPARED  WITH  THEORY  FOR  H/D  =  0.25,  U*.  =  191  FPS 
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Figure  23<c).  VARIATION  OF  EXPERIMENTAL  VALUES  OF  W^/R  WITH  r/D 
COMPARED  WITH  THEORY  FOR  H/D  =  LO  ,  UM  =  191  FPS 


Figurt  23(d)-  VARIATION  OF  EXPERIMENTAL  VALUES  OF  f/y/R  WITH  r/D 
COMPARED  WITH  THEORY  FOR  H/D  =  2.0,  U—  =  191  FPS 
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Figure  23(e).  VARIATION  OF  EXPERIMENTAL  VALUES  OF  y%/R  WITH  r/D 
COMPARED  WITH  THEORY  FOR  H/D  =  4.0 ,  =  191  FPS 


Figure  25.  EXPERIMENTAL  CURVE  OF  DYNAMIC  PRESSURE  IN  MIXING  REGION 

SURROUNDING  THE  POTENTIAL  CORE  OF  THE  FREE  JET,  U*  =  121  FPS 


t/D 


F  lour.  26.  CONSTANT  VELOCITY  CONTOURS  IN  MIXING  REGION 

SURROUNDING  TH  POTENTIAL  CORE  OF  THE  FREE  JET 
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PUNE  STATIC  PRESSURE  DISTRIBUTIONS  •  h/D  =  1.0  TO  2.0 


02 


63-11 


Fijure  28.  COMPARISON  OF  PREDICTED  AND  EXPERIMENTAL  JET- CENTERLINE 
STATIC  PRESSURE  DISTRIBUTIONS  -  H/D  =  1.0  TO  2.0 


-j- — 1_ — | - j - i - i — 4 - j- — ~j — — | — ~4~~ — —J-i- 

i — j — } — t — f — f — I . i — | — \ — S — / 

-4 - i - i - i - j - j- - f — f — -j- . — t - i - 


EiE«EH±E3EEiE^t£2E 


— ; - v - U-- 


^■tTl+Fi/F 


■  ********|*  ■”*****+•*♦•”*  *-*r-***"***V*”  •■•i**”****** 

. [■— — f — — - - - **t- . — — — j- - 

. i . i  i . i  t'""i  n 

- 1 - ; - 1 - 1 - \ - ; - — -J" 


...r - - +-* 

~t - , - r.. 

— r - i - 


3 


*flw***Mt - ♦ . + - 


Figure  31.  VARIATION  OF  GLAUERT’S  PARAMETER  a  vs 

*  =  0.013 
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